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Abstract 

In the context of the Anderson theory of high T c cuprates, 
we develop a BCS theory for Luttinger liquids. If the Lut- 
tinger interaction is much stronger than the BCS potential 
we find that the BCS equation is quite modified compared 
to usual BCS equation for Fermi liquids. In particular T c 
predicted by the BCS equation for Luttinger liquids is quite 
higher than the usual T c for Fermi liquids. 
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1 The anomalous BCS equation 

The equivalent of BCS theory for a Luttinger liquid has 
not formally worked out, despite the relevance of such the- 
ory for the problem of superconductivity in the high-T c 
cuprates, see [1] or the discussion in the following sec- 
tion. Such theory should describe superconductors which 
in their normal state are Luttinger liquids. In this paper 
we develop such theory using constructive quantum field 
theory techniques, applied in many other Luttinger liquid 
problems, see [2], [3]. In particular, we compute in a rigor- 
ous way the BCS self-consistence equation for a spinning 
Luttinger liquid (the Mattis model, see [4]) coupled with 
a BCS anomalous potential. We stress that it should be 
difficult to find similar results using different techniques: 
in fact the Mattis model plus a BCS interaction is not 
exactly solvable, and bosonization or conformal quantum 
field theory cannot be used to compute the correlation 
functions of theories with gap like the present one (see for 
instance [5]). 

We will assumef] that an anomalous self energy has been 
introduced by some external influence in a Luttinger liquid 

1 See [1] pag.209 



where tpu^ a are fermionic creation or annihilation oper- 
ators with momentum^ k = 2m: /L, —L/2 < n < ^j^-, 
quasi-particle index cj = ±1, spin a = ±i; \v(p)\ < e _K ' p '. 
We assume that A > 0. The BCS interaction is described 
as usual, if g > 
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We assume that A e R and A > 0. Note that if A = 
the model is exactly solvable [4] ; the system is a Luttinger 
liquid (it is perhaps the simplest spinning model showing 
Luttinger liquid behaviour) and shows spin-charge sepa- 
ration. 

The ground state energy E (A) depends on A; the BCS 



equation is the extremizing equation 



dE (A) 
dA 



which has 



the form (see [6] for the proof of a similar statement): 
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2 We are considering the Mattis model on a lattice, while the orig- 
inal model introduced and solved in [4] is defined on the continuum 
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If A = the r.h.s of eq.(Q) can be calculated very easily, 
as the hamiltonian can be put in diagonal form by per- 
forming a Bogolubov transformation. In the interacting 
case an exact solution is not available and we compute 
the r.h.s. of eq.(4) writing it as a convergent series for 
A, A suitably small i.e. A, A < e, if e is a suitable num- 
ber << 1 (but no restriction is imposed on their ratio). 
Note that the restriction on the smallness of A, A is done 
only for technical reasons in order to ensure convergence 
of the perturbative series. However it is very likely that 
our bound for the convergence radius are far from to be 
optimal, and instead our results hold also for A = O(l); if 
A = we know that this is true directly from the exact 
solution [4] []. 

Our main result is 

There exists an e > such that, for A, A < e the BCS 
equation (Qj, in the limit j3 —> oo,L — * oo, can be written 

nA A 
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where 77 = /3iA + 0(A 2 ) is a critical index, f = /(A, A), 
l/l < const, and A, a, Pi are positive constants. 

Eq.(^) is difficult to solve in general; we study its non 
trivial solutions for two limiting cases, denoting by cr{x) a 
function bounded by |x| times a constant. 

1) - << 1, which means that the BCS hamiltonian 
is weakly perturbed by the Luttinger interaction. As one 
expects, the gap is given by 



(G) 



so that the gap is exponentially small and the Luttinger 
liquid behaviour of the system plays no role. 

2) j >> 1, which means that the Mattis hamiltonian 
is weakly perturbed by the BCS interaction. We find, re- 
membering that 77 = f3i\ + 0(A 2 ) 



A = J 4[-£-]i[i + a (A) + C r(£)]l 
arj A 
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We see that the Luttinger interaction enhances strongly 
the gap with respect to the A = case if - >> 1 



(8) 



Note also the crucial role of the sign of the Luttinger in- 
teraction. 

2 A model for the high-T c cuprates 

In the previous section we have written a BCS equation 
for a Luttinger liquid and we have seen that, if the Lut- 
tinger interaction is much stronger than the BCS interac- 
tion, there is a large deviation from the Fermi liquid BCS 
equation and in particular the gap is much larger. Can 
this result be applied to the physics of high-T c cuprates? 

The Anderson theory of high-T c cuprates is based on 
the following points^: 

1) the conduction electrons are confined on layers; 

2) thc interaction between electrons in the same layer is 
much stronger than the interlayer interaction; 

3) the normal state of the electrons on the layers is a 
Luttinger liquid; 

4) the interlayer pairing allows Cooper pairs to tunnel- 
ing into an adiacent layer by the Josephson mechanism^ 

The interlayer interaction is assumed to be given by the 
Josephson pair tunneling hamiltoniar^ and some physi- 
cal arguments for motivating this choice and for not con- 
sidering single particle tunneling are given[]. The model 
is then studied by the usual mean field BCS approxima- 
tion^. However by point 3) one has to take into account 
in the resulting BCS equation the Luttinger liquid nature 
of the fermions in their normal state. This is a big prob- 
lem as no theory showing the Luttinger liquid behaviour of 
bidimensional strongly interacting fermions exists^ so that 
the form of a Luttinger liquid BCS equation is essentially 
guessed by replacing in the usual BCS the Fermi propa- 
gator with the one dimensional Luttinger model propaga- 

As the above theory is a mean field theory, its results 
should be, as usual, independent from the dimensions; this 



4 See for instance pag. 46-55 of [1] 
5 See also pag. 308 of [1] 



6 See for instance eq. (7) pag 321 of [1], or eq. (9) (in d = 2) below 
7 See for instance the considerations based on the holon-spinon 

coherence at pag. 50 of [1] or in the Clarke reprint [o] 

8 see eq.(16) pag. 218 of [1] or the considerations after eq.(10) 

below 

9 In [1] this is essentially derived from experimental results 
10 See pag.213 eq.(12) of [1] 



3 The role of the lattice should play no role. 
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means that we can expect that coupling d = 1 or d = 2 
Luttinger liquids and performing a BCS approximation on 
the pair tunnelling interaction, the results are qualitatively 
the same. In other words we expect that the predictions of 
a BCS theory for a superconductor whose normal state is 
a Luttinger liquid are qualitatively the same whatever its 
dimensions are. This is just what happens for the usual 
superconductors whose normal state is a Fermi liquid; the 
BCS approximation leads for instance to a gap or a critical 
temperature exponentially small in g in any dimension. 
The crucial and very non trivial assumption is that there 
exist d = 2 fermions with a Luttinger liquid behaviour in 
their normal state; but, if they exists, the results of a BCS 
theory of coupled Luttinger liquids should be essentially 
indipendent on their dimensions. 

Considering one dimensional coupled Luttinger liquids, 
the BCS equation can computed in a rigorous way. The 
Luttinger liquids can be described by the Mattis model 
hamiltonian eq.(2), giving an extra chain index to the 
Fermi operators i = a,b. The two Mattis hamiltonians 
will be called H a and H^. The pair-hopping hamiltonian 
Hint is, following 

ff « = -«[iEWWiJ (9) 
[yE^,- Ur i,A, u ,i,J +cx ' 

The total hamiltonian is H a + H\, + Hi nt . Let we write in 
eq.® 

— V ib + ib + 

as 

iE<l^Wi, >+ (10) 

k.LO 
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The BCS approximation consists in replacing Hi nt with 
Hint obtained neglecting the terms bilinear in the "fluctu- 
ations" i.e. in the second addend in eq. ©. We obtainQ 

k,u 

n See eq.(7) pag.320 
12 See eq.(16) pag.218 of [1] 
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and self-consistency requires eq.(^). Replacing Hi nt with 
^int P nas * ne en?ec t that the model is described by two 
indipendent hamiltonians, H = H a + Hi,, each one given 
by eq.(l). 

The BCS equation for coupled Luttinger liquids is then 
given by eq.(5). In the range of parameters physically 
reasonable (see point 2) above) i.e. if A >> g and noting 
that, as usual, T c is proportional to A (see below) we find 

~ — f - I 1 3 1 1 

that T c ~ Be ° 8i "i" i.e. much higher than T c for 

normal superconductor T c ~ B'e~s. 

In order to compare our result with [1], note that the 
r.h.s. of eq.(4) computed^] in [1] is similar to our eq.(p2|) 
below (see especially the presence of the wave function 
renormalization) but ah ~ A-y-(@i > '+--) h is replaced by A 
i.e. the anomalous flow of the BCS gap is neglected. Such 
effect is on the other hand crucial for our analysis: all the 
conclusion drawn from eq.(5) in the preceding section are 
based on the fact that j3\ ^ 0. In fact it is this anomalous 
enlarging of the gap due to the Luttinger interaction which 
produces a much larger solution of the BCS equation. On 
the other hand the gap renormalization in Luttinger liq- 
uids is a well known fact, see for instance the gap of the 
XYZ model, [9]. 

Finally we stress that one can try to study directly a 
model of coupled chains, see [7], without any approxima- 
tion, and it could happen that the model is not really 
described by our BCS approximation. This has not real 
relevance for our analysis, as our aim was just to find a 
BCS equation for coupled system of fermions on planes 
with a Luttinger liquid behaviour, and we make a BCS 
computation in d — 1 using the fact that a mean field the- 
ory like BCS should be insensitive to the dimensions. So 
it could be that our BCS equation eq.(5) could be applied 
to coupled planes and not chains, and on the other hand a 
detailed analysis of the model given of eq.(10) could give 
no insight on the problem of coupled planes as the be- 
havior of the system depend on the dimensions. This is 
in fact just what happens in the usual BCS theory: the 

13 The above computation follows the Anderson gap derivation, 
see chap. 7 in [1], but for semplicity we have neglected the electron 
phonon interaction, see eq. (21) pag. 219 of [1] 

14 see eq.(12) pag.213 
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BCS equation is qualitatively the same in any dimension, 
but in d = 1 the mean field approximation is not correct. 
The real question is if really bidimensional Luttinger ex- 
ists, but this question is not addressed here and we refer 
to [1]. 

3 Renormalization group analysis 

We discuss now how to compute the r.h.s. of eq.([|). We 
introduce as usual a set of Grassman variables ihi 

k = (k ,k) and k = k = 27r ("°+ 2 " 1 ) ; if no ,m are 

integers and -| < m < (L - l)/2, -§ < n 2 < 
if 1//3 is the temperature. The Grassmanian integration 
Pldip) is defined by the anticommutative Wick rule with 
propagator 

g u (k) = — 

—iko + ujk — pp 

In general we denote by J {T>ifie~ ^ h ^ ^ £ } the 
Grassmanian integration with propagator h(k), where 
J dk — Sfci m particular 

P(#) = W n e - $ dU i^ { - iko+ " k - pr) ^} 

Then we write the r.h.s. side of eq.(d) as a functional 
integral 



1 



i 



(ii) 



where 7(V>) = AF + AP and 



4 4 
i— 1 i— 1 

E: i^T i/iIT :: i^T r^IT : 



We evaluate the above Grassman integral using Wilsonian 
renormalization group techniques. It is convenient to write 
k = k' + lopf, where k' is the momentum measured from 
the Fermi surface. 



We decompose the integration P(dip) into a product of 
independent integrations. This can be done writing 

— oo 

g w {k' + uj Pf ) = g%(& + u PF ) + gift + uj Pf ) (12) 

h=0 

with 



gt( k ' + uj pf) 



fh(k') 



-iko + Lok' 

and/i(fc') = 1— x(fc')i x(k') = x(l^'l) is a smooth compact 
support function such that = 1 f° r W\ 7 _1 an d 

x(k') = for > 1, if 7 > 1; moreover for ft, < 

h{k')=x(l- h k')~x{l- h+1 k') 

is a smooth compact support function non vanishing only 
for 7^ _2 < \k'\ < ry h . Then g^(k' + u>pf) is the ultraviolet 
part of the propagator, while X)/c^o5w(^') * s the infrared 
part. Note that, from the compact support properties of 
g^,, the sum in eq.(12) is from to hp, where ^ hli = 5, as 

|fco|>|. Let be C^" 1 = Eth, /*• 

The ultraviolet integration is somehow special (and es- 
sentially trivial for the presence of the lattice) and we will 
not discuss it here, see [2]. If A = the infrared integra- 
tion can be done by performing the well known Bogolubov 
transformation to diagonalize the BCS hamiltonian. If 
A 7^ the BCS gap and the wave function renormalization 
have a non trivial RG flow so that we have to perform 
a different Bogolubov transformation at each iteration of 
the RG. We set Z = 1: once the fields ^ (0) , ■ ■ • , i> [h+1) 
have been integrated we have: 



Uf=±l 

v' l (v / z^V ( -' l) ) 



- f dk'c h z h ^ h)+ g^ (kT 1 ^ 

<J k'io k'u 



} 



(13) 



if 



^ = (4^ j 

k ,UJ k -\-LOpF ,UJ, 2 — K, — LUpF , — UJ,— 2 

and g( /l )(fc') _1 is defined by 

(—iko + ujk' o~h(k') \ 
<Jh{k') ~ik a - ujk' J 

V h is called the effective potential at scale h and is given 



by 



(<h)\ 



oo p n 
n=2 J i=l 



(14) 
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wZ(k 1: ...X)f\^ h)ei {(y^+w)) where Q{h x -f } x - is , defined as i n eq -@' with ^~ 1 

nv i ^H^,^ istcadof h, V h = £V h + {1 - £)V h , 



i=l i=l 



We define a localization operator £ extracting the rele- CV h = j h nh F h + (a^ — z-^)F h + g 2 ^-Ft 1 + 34 /i^i 1 
vant or marginal part of the effective potential V h : 

i) If n > 4 then £W* = 0; and 

ii) Let be n = 4. In this case £^4 = unless Zh-i{k') = Z h + C h ~ 1 Z h Zh 



^2^ =1 6iLUiPF — 0, X)i=i e » = in which case the action 
is non trivial and it is given by 



^_i(fe')^-i(A;') = Z h a h {k') + Z^C- 1 ^ 
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This means that we extract from the effective potential the 
£W A (k 1 +u>ipp, ....) = W 4 (ujiPf, ■■■) (15) terms leading to a mass and wave function renormaliza- 
) if n — 2 then if y"V e • — tion. Now one can perform the integration respect to ip( h \ 

rescaling the effective potential V h (ip) — ^ h {\p^~^) 
CW^ik't +u;ip F ,k' 2 + uo 2 Pf) = \W 2 {u)ip F ,uj 2 pf) + and 1 



uj\E{k' + LJipp)dkW 2 (ljiPf, uj 2 pf) + k° 

d ko W^U lP F,LU2PF)] (16) 



with v h = #>-n h) 4 = -^-(a h -z h ), A ilh = (-J^-) 2 gi,h 



while if ^ then J 1 _ z h _i " ft » - v-" «i,ft - vz k _, 

and Vk = Vh, $h,^h- We can rewrite eq.(18) as 



CW 2 (k' 1 +ujipF,k 2 +lu2Pf) = W 2 {ujiPf,uj2Pf) 



We can write then the relevant part of the effective 
potential as: [vijj {h) e~ I dk'^+ f- 1 z h ^ 1 (k')g^ h - 1 Hk'r 1 ^ k f : y 



CV h = j h n h F!; + Sh F2 + z h F£ + a h FZ + g2, h F> l + g 4 , h F> 1 -v\^r^^ 

(17) 6 1 

wnere and the integration of ip^ has propagator 



F n = Y'[dk'M ( i h)+ ^ h) - , 1 r 

f£ = Y I dk'^ h)+ ^ 
F 2 h = [[f[dh6(J2ki) E 

i— 1 i uj,<j,<t' 



dk'e^'^flik^g^ik')^ 



(<h) + 

- -fee 
+W p J? ,u;,i r -fc'- W p f .,-a>-A •• with ^^(A;') given by 



u,a,a< A h \ -ah-i{k') (-iko-uik')) J 

r^(p),+ ^(p),- ^(p),- j 

f 4 = /in &) s and = z >'-i(°) and 

z— 1 i uj,o,(t' 

[v ,<p),+ ^(<M,- ][v ,^),+ ^(p),- ] 7 z 1 r__£_____ki 

where i — v, a, Q, f v — 1, f a — u>k , Jq — *^o- More rp^ e regu ^ Q £ integration is in the same form as eq.(p"3|) 

0VCr 9W = ^°) A + °( A2 )' ^ CA2 ' So = A + °^ AA )' with h replaced by ft - 1, and we can iterate. 
GO) z o 0(A), no O(AA). We write eq.(|13|) as. wg explain the main motivations of the integration 

jy^p(<h) g - J dfc'^7 h)+ c h z h _ 1 (fe')e <h_1) (fc') _ Vfc7 h)_ procedure discussed above. In a renormalization group 

approach one has to identify the relevant, marginal and 



(18) irrelevant interactions. By a power counting argument 
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one sees that the terms bilinear in the fields are relevant 
and the quartic terms (or the bilinear ones with a deriva- 
tive respect to some coordinate acting on the fields) are 
marginal. However there are too many kinds of marginal 
terms, depending on the labels uji and on each fields, 
so that their renormalization group flow seems impossible 
to study. However (see [2], [3] for a similar procedure) 
the power counting can be improved and many marginal 
terms are indeed irrelevant; in particular all the marginal 
terms with four or two fields with ^ i &j ^ are indeed 
irrelevant. The reason is that such terms are generated 
contracting at least a non diagonal propagator and such 
propagators are smaller than the diagonal ones by a fac- 
tor (Jh7~ h , see eq.([l9]) below; this will be sufficient for 
improving the power counting, (see the last paper in [3] 
for the proof of a similar statement in the XYZ chain). 
Moreover also the marginal terms with J^i £i^iPF ^ are 
irrelevant, by momentum conservation considerations. In 
fact if Yli £ i u, iPF then the momenta of the fermions 
cannot be all close to the Fermi surface; mathematically 
this means that, for the compact support properties of the 
propagators, there is an h such that all scattering process 
involving fermions such that J2i £ i UJ iPF ^ with scale 
h < h are vanishing. 

The relevant terms are of two kinds; the v terms, reflect- 
ing the renormalization of the Fermi momentum, and the 
a terms, related to the presence of a gap in the spectrum. 
The presence of the v terms is due to the renormalization 
of the chemical potential, and in general one introduces a 
counterterm in the hamiltonian to fix the Fermi momen- 
tum, sec [2]. In this case however there is no necessity of 
adding this counterterm; roughly speaking, fi can vary in 
the gap whithout changing the Fermi momentum i.e. the 
position of the singularity of the propagator. This is a cru- 
cial point: if one had to put a IS.- dependent counterterm 
in the hamiltonian, then considering dE Q^ one would be 
forced to derive also such counterterm, and a much more 
complex BCS equation would appear. Regarding the other 
relevant term, they are related to the BCS gap genera- 
tion. However due to the interaction the BCS gap has a 
non trivial flow, so that one has to perform different Bo- 
golubov transformations at each integration. 

Regarding the marginal terms, there is an anomalous 
wave function renormalization which one has to take into 



account, what is expected as if A — the theory is a 
Luttinger liquid. In general the flow of the marginal terms 
can be controlled using some cancellations due to the fact 
that the Beta function is "close" (for small u) to the Mattis 
model Beta function. In eq. ( pfj| ) we write the propagator as 
the Mattis model propagator plus a remainder, so that the 
Beta function is equal to the Mattis model Beta function 
plus a "remainder" which is small if o~hJ~ h is small. 

Let be h* = inf h {7' 1 > \<r h \}. Note that, if h* is finite 
uniformly in L,[3 so that \o-h*-i\l~ h +1 > 1 one has 



\g <h '\x)\< 



i 



C M l h * 



Z h . \ + {^\x\) M 
Moreover if h > h* we have 



Iff" (*)l < 



1 



C M l h 



Z h l + (j h \x\) 



M 



and 



\JUJ.—UJ 



C M l h 



(19) 



Z h l h 1 + (7 ft |^|) M 
Moreover for h > h* the bound for the non diagonal 
propagator has a factor more with respect to the di- 
agonal propagator. This is the reason for which the quartic 
terms with J^. e% ^ are irrelevant, despite dimensionally 
marginal. Finally 



(x-y) = gi L (x-y) + cL(x-y) (20) 



with 



-v)= dk' 



Ah' x 



h(k') 



Zh —iko + wfc' 
which is just the propagator "at scale h" of the Mattis 
model, and the other term verify the bound of 9u U (S;y) 

with an extra factor ^V-. 

7 

We see from the above bounds that the propagator of 
the integration of all the scale between h* and /tgPlhas the 
same bound as the propagator of the integration of a sin- 
gle scale greater than h* ; this will be used to perform the 
integration of all the scales < ft* in a single step, i.e. inte- 
grating directly ip( <h '. In fact 7^ is a momentum scale 
and, roughly speaking, for momenta bigger than j h the 
theory is "essentially" a massless theory (up to 0(ahj~ h ) 
terms) while for momenta smaller than j h is a " massive" 
theory with mass 0(7 ). 

15 Of course if h* < hp there is no such integration. 
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One can prove that the effective potentials V h are 
well defined, if the running coupling constants are small 
enough. More precisely, let we write eq.(14) in coordinate 
space and let be W% the corresponding kernel; it holds 
that 

Lemma: Assume that h* is is finite uniformly in L,[3 
and that for any h > k > h* there exists an e such that 
|4| < e and |2£±l| < Y" e , < Y"^ with c a ,c b 

postive constants. Then there exist a constant C such that 



\W*\\ < Nf3eC~f- k{ %- 



■2) 



The proof of the above lemma is an immediate modifi- 
cation of ones existing in literature, see in particular [3]. 

In order to prove that the effective potentials are well 
defined we have to show that the above conditions of small- 
ness in the above lemma on the running coupling constants 
are verified. 

The beta function for Vh is 

Vh-\ = lVh+ Ph+Ph 

where f3h is the contribution to Vh obtained setting uy = 
0, < k' < h, so it is exactly equal to by the parity 
properties of the Mattis model, and 

\Ph\ < C^X 2 

is the remaining part. Iterating the above relation we find 
= 1 W 7 " l * l A 2 ]T K| < C|A| 

k=0 k 

as a k ~ A^ lk , m = -/3iA + 0(A 2 ), see below. 

The Beta function can be written, for > h > h* . : 

M,h-i = A/j + G\ h + Gl' h 
M,h-i = A/i + G\' h + Gl' h 
o~h-\ = &h + Gl' h 
5h-! = S h + G]' h + G 2 5 ' h 



z h 



1 + Gl Jl + G 2 / 



(21) 



where (i — 2, 4) 

a) G\' , G\ h and G\ h depend only on A i)0 , 6 ; ...Xi,h, Sh 
and are given by series of terms involving only the Mattis 



model part of the propagator g k L (x — y), so they coincide 
with the Mattis model Beta function 

b) G 2 a ' h , G 2 ' h , G 2 ' h , G 2 z ' h are given by a series of terms 
involving at least a propagator C\ u (x — y) or g k _ u (x — y) 
with k > h. 

By a simple explicit computation 

Gl< h = \ 2 [fo + G h z } 

G 1 * = \ h v h \/3 1 + G* r ] 

with (3 1 ,fa>0 and G h z , G h a = 0(X h ). Moreover G]' h , G 1 / 1 
coincide by definition with the Mattis model Beta function, 
and it was proved in [2], [3] that it is vanishing at any 
order, i.e. 

G i,h = G i,h = 

Finally as \G 2 ' h \, \G 2 ' h \, \G 2 z < h \ < Ke 2 \a h \ 1 - h , one finds, 
for h > h*, 

Wh-i - Aj,o| < ciA 2 \5h-i - 5 \ < ciA 2 

and a h ~ (&h Vlh , Z h ~ for h > h* , with 

?7i = -/3iA + 0(A 2 ), 772 = /? 2 A 2 + 0(A 3 ). As usual in 
models to which the RG is succesfully applied the flow is 
essentially described by the second order truncation of the 
beta function. This shows that it is possible to choose A so 
small that the conditions of the above lemma are fulfilled. 
From the definition of h* and the fact that ah ^ (A)"/ 11111 
it follows a h * = A 1-/3iA+ °( a2 ). 

As we said the integrations of the ip( <h ' (if h* > hp) 
is essentially equivalent to the integration of a single scale 
h > h* , so it is well defined by the preceding arguments. 
If h* < hp there is no such integration, and the last scale 
to be integrated is hp; from this consideration one obtains 
easily that the critical temperature is proportional as usual 
to the gap amplitude. 

An expansion for the two points Schwinger function can 
be derived in a standard way [2] and from the proof of the 
convergence of the expansion for the effective potential 
one obtains easily the convergence of the series for the 
Schwinger function. We can write the r.h.s. of eq.(4) as 



E 



Oh 



/i-max[/i* 



(27T) 2 

L(3 4^ Z h k 2 + k /2 +a 2 h 



E 



h(k') 



[l + \S h (k')] 



(22) 
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where Sh(k') is a convergent series bounded by a constant; [9 ] A.Luther, Phys. Rev. B, 14,5, 2153-2159 (1976) 
from the above expression one can easily derive the BCS 
equation eq.(5) as well as the critical temperature. 

Note the crucial role of the renormalization of the BCS 
gap ah) it is sensitive to the sign of the interaction and 
it eliminates or enlarges the singularity of the r.h.s. of 
eq.(5); neglecting such renormalization one obtains com- 
pletely different results. In fact our model belongs to the 
class of universality of the massive Luttinger model, for 
which it is well known that the bare mass A is renormal- 
ized by the interaction to be given by A 1 ~ n , n = O(A); 
other models belonging to this class are the XY Z chain 
or the Yukawa2 model, see [3] , [9] . Note also that no role is 
plaid in the above analysis by the spin-charge separation; 
in fact \v c — v s \ = 0(A) and such effect is incorporated in 
the term S h in eq.(p2|) 
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